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ABSTRACT 

We have used RR Lyrae and Blue HB stars as tracers of 
the old Galactic halo, in order to study the halo structure 
and the galactic rotation as a function of height above the 
plane. Our sample includes 40 RR Lyrae and 80 BHB 
stars that are about 2 to 15 kpc above the plane, in a 
roughly 250 deg^ area around the North Galactic Pole 
(NGP). We use proper motions (derived from the GSC-II 
database) and radial velocities to determine the rotation 
of the halo. From the whole sample the motion appears 
to be significantly more retrograde than the samples in 
the solar neighborhood, confirming Maiewski ( 199% re- 
sults and our own preliminary results bas ed on 1/3 th e 
present sample (iKinman et al. 2003; Scagn a^t al.l2003l) . 
However, the better statistics has now revealed the likely 
existence of two components, whose characteristics need 
an accurate analysis of systematic errors on the proper 
motions in order to be assessed in detail. 

Key words: Kinematics, Galactic structure, RR Lyrae, 
Blue HB stars. 



1. INTRODUCTION 

The kinematics of halo stars in the solar neighborhood 
(Z < 2 kpc) is well established: mean heliocentric 
rotation velocity (V) ^ —220 ± 10 kms~^ (from 
RR Lyrae an d K gi ant stars, cf. 'Martin & Morrison' 
[1998; Lavde netalJ IT9961 [Chiba & Yoshii 1998; 
pambis & Rastorguev 2001). 

From halo stars observed at larger distances {Z > 5 kpc), 
the kinematics of the high halo is less clear: the mean ro- 
tation was found retrograde {(V) = —275 ± 16 km s~^) 
by Maiewski ( 1993) and 'Maiewski. Munn & HawlevI 
(Il996i) in SA 57 near the NGP, and confirmed by 
iGilmore. Wvse ifeNorris (.2002) in other directions. 
Howe ver, othe r estimat es calculated from o rbital param- 
eters ("Carnev 1999, iChiba & BeersI l2000h . as well as 
fSirko et al. (2004) recent analysis of very distant BHB 
stars from the SDSS, do not support any significant 
halo rotation. The present work attempts to resolve this 
discrepancy; it grew from earlier studies of halo stars 
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in the North Galactic Cap (iKinman et alJll99 6^ which 
confirmed the streaming motion (in the W vector) that 
Majewski found for his subdwarf sample in SA 57. 
So the questions are: /) Is the high halo in retrograde 
rotation? ii) Is there any halo substructure that may 
mimic retrograde rotation in some directions? 



2. OUR DATA 

Our sample of halo tracers consists of 80 blue HB (BHB) 
and 40 RR Lyrae (RRL) stars spread over an area of 
approximately 22° x 12° around the NGP, at a distance 
Z Cf. 1.5 to ^ 16 kpc above the galactic plane. This 
area is the combination of 6 POSS fields where proper 
motions from the GSC-II database were measured (see 
later), and is shown in Fig. Q We refer the reader to 
Kinman et al. ( 2003) for details on the target selection. 
For these stars we have: 

• Absolute Magnitudes (±0.2 mag) =J> dis- 
tances. In particular, E(B — V) values are from 
ISchlegel. Finkbeiner & DaviesI (Il998h : for the BHB 

stars. My have been estimated from the Mv vs. {B — V) 
relation given by Preston, Shectman & Beers (1991) 
adjusted so that My = 0.6 at {B - V) = 0.2; for the 
RRab stars, we used My = 0.22 [Fe/H] + 0.93 if 
\Fe/H] was av ailable, or My = -1.619 logP + 0.20 
(iKinmanl l2002h if only the period was available, or 
else My = 0.6. For the RRc stars it is assumed that 

My = 0.6. 

• Proper Motions (formal r.m.s. errors < 3 masyr"'^). 
These are based on the plate material used for the 
construction of the GSC-II cata logue using the method 
described i n ISnagnaet 311(119961) . 

• Radial velocities (±10 — 50 kms^^), obtained at the 
4m-RC (KPNO) and 3.5m-LRS (TNG). 

From the above data we have derived heliocentric space 
ve locity components U, V, and W using the program 
by iJohnson & SoderblomI (Il987l updated version for the 
J2000 reference frame), with a further update of the trans- 
formation matrix derived from the Volume 1 of the Hip- 
parcos data catalogue. Finally, the heliocentric UVW ve- 
locities have been corrected adopting the solar motion 
(U,V,W)0 = (10.0, 5.25, 7.17) km s-i with respect to 
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Figure 1. The P OSS fields at the NGP where we have 
performed the present analysis. Right ascension, along 
the X axis, is expressed in arc degrees. 

the LSR, from lDehnen & BinnevHl998h . 



3. ERRORS 

An error in Aly will not have a large effect on the derived 
U and V velocities (a total uncertainty of as much as 
±0.2 mag would correspond to ±10% in the distance and 
hence V). The U and V velocities are nearly independent 
of the radial velocity in the direction of the NGP. 
In this case, the greatest uncertainty in V will come from 
errors in the proper motions. 

Therefore it is important to test the systematic errors 
in the proper motion (pm) data using the QSOs and 
galaxies present in the observed fields. Proper motions 
have been derived for 65 QSOs and 50 galaxies, omitting 
a few objects whose pm or errors exceed 10 masyr"^ 
(iKinman et alJl2003h . These 115 objects, which should 
have zero proper motion, have the following average 
GSC-II pm: 



fia = -0.109 ± 0.122 masyr-1 
^ls = -0.422 ± 0.162 masyr-i 



(1) 



Presumably these systematic errors (< 1 mas yr^^, com- 
parable to those of the pm measured by Hipparcos for 
much brighter stars) apply to our program stars over the 
same sky area and a similar magnitude range. 
A preliminary check shows systematic errors of similar 
size in our individual fields. Although these individual 
field errors are necessarily more poorly determined, we 
have used them to correct our BHB and RRL motions in 
this preliminary analysis. The distribution of V shows lit- 
tle systematic trend with either position or distance (see 
Fig. 13, so these corrections introduce little or no bias. 
Currently, the accuracy of the GSC-II proper motions rep- 
resents the best available for the large-field surveys that 
are needed for studies of Galactic structure; higher preci- 
sion is necessary but must await dedicated space missions 
such as GAIA. 



4. RESULTS 

In calculating the UVW vectors we put the radial veloc- 
ity equal to zero (with an error of ±150 kms^^) if no 
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Figure 2. Rotational velocity V after correction for pm 
systematic errors in each field, as a function of position 
and distance. Different symbols represent different V in- 
tervals, and the size of the symbol is proportional to the 
parallax. 



radial velocity was available. In such cases (7 RRL stars 
at Z > 4 kpc), the U and V vectors should be scarcely 
affected but the W vector is discarded. For a better char- 
acterization of our stars, we have also trimmed 10% of 
the most deviating stars from our sample when estimating 
mean values or distributions. We have found that trim- 
ming has little effect on the mean values of U, V and W, 
although it does reduce the velocity dispersions (Ttj, cry 
and iTw In Tabled we compare the mean heliocentric 
UVW values of our entire sample and of the subsample 
at Z > 4 kpc with those found by Martin & Morrisoi^ 
(|1998^ for their Halo2 sample of local RRL stars (in all 
cases 10% of outliers were trimmed). We also show the 
W-U, W-V and U -V plots in Fig. E] 
A KMM test (cf. Ashman, Bird & Zepf 1994") on the en- 
tire trimmed sample gives about 90% probability that the 
sample is formed by two groups, one containing ^ 56% 
of the stars with estimated {V) — —316 ± 8 kms~^, and 
the other containing ~ 44% of the stars with estimated 
{V) — — 177±9A:mi"^, same dispersion cr = 61 kms^^. 
We show in Fig. 0]the distributions of V as a function of 
the distance Z above the Galactic plane, where the bi- 
modal shape is evident at all distance intervals as well as 
for the entire sample. 



5. CONCLUSIONS 

Our present results, based on a 3 times l arger sample 
than our previou s preliminary analysis ( Ki nman et alJ 
I2003L ISpagna et al..,.2003.) . seem to lead to the following 
conclusions: 

• For Z < A kpc, the mean rotation V is close to 
that found for the halo in the solar neighborhood (see 
corresponding panel in Fig.0}. 

• For Z > 4 kpc, the mean rotation is significantly 
more retrograde - in agreement with Majewski's earlier 
finding. This seems to be due to a group of very retro- 
grade stars, whose relative importance increases with 
Z with respect to the dissipative halo component, and 
is particularly significant in the range Z = 4 — 10 kpc 
(see Fig. |4}. Therefore, this is likely to be a local 
substructure, maybe associated with an accretion event. 
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Table 1. Space motion vectors for the present halo stellar sample, and for a local halo sample for comparison. 



U V W (Tu (Tv iTw No. of Stars Source 



-12±15 -256±9 -10±9 153 92 91 108 (101 for W) Present sample 

+4±15 -264±9 -20±9 158 91 95 82 (75 for W) Present sample at Z > 4 kpc 

-1±26 -219±24 -5±10 193 91 96 84 Martin & Morrison (1998) local sample 
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Figure 4. Histograms of rotational velocity V as a func- 
tion of distance Z above the Galactic plane: all distribu- 
tions suggest a bimodal structure. 
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since ISirko et aP (12004') analysis of a large number of 
distant BHB stars spread over the sky shows little devia- 
tion from the solar neighborhood (dissipative halo) value. 

It will be interesting to extend the work to halo stars in 
Anticentre fields. Possibly this will allow us to detect 
gradients in the V motion and discover whether this 
is a local effect limited to the NGP or part of a larger 
systematic effect. 

The errors in the GSC-II are small enough to make our 
results worth consideration. On the other hand, they are 
large enough to make it clear that GAIA proper motions 
are essential for a detailed knowledge of the Galactic 
structure. 



Figure 3. Plots W-U, W-V and U-V of the entire stellar 
sample, including the "outliers " that were trimmed out 
in the previous considerations. Triangles and circles rep- 
resent RR Lyrae and BHB stars, respectively. The mean 
heliocentric rotation velocity (V) ~ — 220±10^ms~^ is 
reported in the V plots (shaded line). Note that in the dif- 
ferent velocity planes the contamination from disk stars 
is always quite small. 
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